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We compute derivatives of thermodynamic quantities with respect to /i, at = 
for 2 and 3 flavors of degenerate quark masses. This allows us to estimate the 
phase transition line in the T, ^ plane and quantify the influence of a non vanishing 
chemical potential on the equation of state by computing lines of constant energy, 
pressure and density. Moreover we evaluate the order of the QCD phase transition 
by measuring the Binder Cumulant of the chiral condensate. This gives access to 
the chiral critical point on the phase transition line. 



1 Introduction 

To understand recent heavy-ion collisions, it is mandatory to study the QCD 
phase diagram for high temperatures and small chemical potentials. It is 
well known that the order of the phase transition is strongly dependent on 
quark masses and chemical potential. In the space of two degenerate light 
quark masses rriu^di one heavier mass and the quark chemical potential 
/i, one expects a critical surface, which bends over the quark mass plane and 
separates the regime of first order phase transitions from the crossover regime 
as shown in fig. 1. 

To determine the order of the phase transition we compute the Binder 
Cumulant of the chiral condensate, B4, at /i = for several quark masses at 
the pseudo-critical coupling j3c{m). We define /3c as the peak position of the 
chiral susceptibility. is given by 

This quantity is a renormalization group invariant quantity, with a universal 
value at m^mcrit- The universality class is that of the 3d Ising model with 
^4 = 1.604EJi3'lI. The critical surface is then given by the surface of constant 
B4 = 1.604. 

To explore the regime of /i 7^ we compute derivatives of transition 
temperature, pressure, energy density, quark number density and Binder Cu- 
mulant with respect to fi. Direct Monte Carlo simulations for /i 7^ are not 
possible due to the sign problem. We therefore use a reweighting method at 
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Figure 1. Phase boundary for the first order regime of the thermal phase transition for 
non- vanishing fj, (left), and the line of second order phase transitions for standard staggered 
fermions at = (right). 



fj, = based on Taylor expansion of the fermion determinant and all observ- 
ables up to order /i^. To be more precise, the expansion is given in terms of 
fi/T as the lattice chemical potential is given in units of the lattice cut-off, 



2 Reweighting in quark mass and chemical potential 



Ferrenberg and Swendsen's reweighting method is a very useful technique to 
investigate critical phenomena. Multi parameter reweighting was first applied 
to the problem of finite density QCD in ref. 4. We ujif the Taylor expanded 
reweighting formula up to order /.t^, which is given by 



((Oo + Oi^i + 02^i'^) exp {Ri^L + R2H^} exp {-A5g}) 
{exp{Rifi + R2fJ.'^}exp{-ASg}) 



Here we have 



Ri — 



Nf d'liidctMifj.) 



Ail 



and 



M=0 



(2) 



(3) 



Ai=0 



We calculate the reweighting operators Ri using stochastic estimators. One 
can easily deduce the corresponding reweighting operators for mass reweight- 
ing. We have control over the sign problem through the odd reweighting oper- 
ators R2j+i, which are purely imaginary. The phase, Q, of the fermion deter- 
minant, det Af = I detM| exp{i0}, is in leading order given by O = filmRi. 
From this we find increasing phase fluctuations for increasing volume, number 
of flavors, n/, and for decreasing m and T. 
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Figure 2. The transition line in 2-flavor QCD together with hnes of constant energy, nuclear 
matter density and the freeze-out temperature (left) and a sketch of the phase diagram in 
3-fiavor QCD (right). 



3 The transition to the QGP and the equation of state 

In 2-flavor QCD the transition from a hadron gas to the QGP will be a 
continuous but rapid crossover for non vanishing quark masses and small 
values of /i. In order to estimate the_transition line we performed simulations 
with 2 flavors of improved fermionsO on a 16^ x 4 lattice and quark masses 
of am — 0.1 and am — 0.2. From the peak positions of the Polyakov loop 
susceptibility and the chiral susceptibility we find cPf3c/dfJ.^ — —1.07(24) and 
(Pfic/diJ? — —1.10(68) respectively, which agrees with results of ref. 7. To set 
the physical scale we use the string tension and obtain a{d(3/da) = —2.08(43) 
at {f3,m,) — (3.65,0.1). Combining these results yields 

(fTc _ 1 (f(3c (dp" 

Derivatives of the pressure p and the interaction measure e — 3p are related to 
the quark number density Uq and quark number susceptibility Xs = duq/djiq 
in the following wayu: 



Tc-j^-~ir.-l^[--] --0.14(7). (4) 



9(p/r4) _ 1 d\nZ Uq 

duq ~ VT^ dfiq ~ T4 



(5) 
(6) 
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We obtain T^d'^{p/T'^)/dfil « 0.69 and T'^d'^{e/T^)/d^il « 10.6 at /3c and 
am = 0.1. In the RHIC regime oi (iq/T^ « 0.1 the deviation ofp/T'^ and e/T* 
from results at /i = thus is only a 1% effect. From the second derivatives of 
p and e with respect to /i together with the derivatives with respect to T, we 
calculate the lines of constant pressure and energy density, which are within 
errors parallel to the phase transition line: 



The quark number density is zero at fi — 0. For /i 7^ It can be estimated via 
n,a3 = fiqaxsa"^- This is njT^ = 0.693(5)/i,/T and n^/T^ = 0.490(4)^5/r 
at am = 0.1 and am — 0.2, which translates into roughly 9% and 6% of nuclear 
matter density at the RHIC point. All these results are sunpaarised in fig. 2 
(left). In addition we also show the freeze-out temperature^. This suggest 
that at sufficiently large /is there is the chance to observe experimentally a 
strongly interacting hadron gas phase, whereas for small hadrons seem to 
freeze out right after the phase transition from QGP to the hadronic phase. 

4 The chiral critical point 

To determine the critical surface shown in fig. l(lcft) we performed calcula- 
tions with much lighter quark masses. For 3 fiavors of improved fcrmionaS, a 
mass value of am = 0.005 and a volume of 12'^ x 4, we measured the reweight- 
ing operators R2 needed for reweighting in m and fi. At present we have 
five /3-values, with a total number of 6100 trajectories. For am — 0.01 and 
V = 16^ X 4 we use Ri = ipip for mass reweighting. To determine the critical 
mass value m, we compute B4 as a function of m. The two different volumes 
have an intersection point near the value of B4 = 1.604, i.e. 3d Ising universal- 
ity class. Due to the large errors we give an upper bound for the critical mass 
only, which is am < 0.0075, or in terms of the pion mass fhps < 190MeV. 

For increasing quark chemical potential fiq we find a decreasing Binder Cu- 
mulant. From the two partial derivatives dB^^/ d{am) and dB4/d{a^fi^), which 
we get from straight line fits of the reweighted data, and the assumption that 
dB^/d^am) is constant in am, one can compute the quantity d{am)/d{a^ ^^). 
The first derivatives dB^/ d{aii) and d(am) / d(a^) vanish because of symmetry 
reasons. A jackknife analysis yields a^^dfh/ d{ii^) — 0.82(23), or equivalently 
Tdm/d{fi^) = 0.21(6). We thus find that the critical quark mass, fh, increases 
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0.07(3). (8) 



with increasing chemical potential as anticipated in fig. 1 (left), 

fh{fi) « TO3 + 0.21 n^/T. (9) 

Here 7713 is the critical quark mass for 3 degenerate flavor at /i = 0, which 
turns out to be roughly twice as large as the physical light quark masses, 
7713 « 2mFj^''J^ . This result is shown in fig. 2 (right). 

Based on these results we can give a first estimate for the location of 
the chiral critical point, i.e. the second order endpoint in the phase diagram 
shown in fig. 2 (left). To do so we use 'mF^^'' /T = 0.016, and assume that the 
curvature of the critical surface is the same for 3 and 2+1-flavor. In eq. 9 we 
should then replace by m2i ~ ^'^^ud''^ which can be deduced from fig. 1 
(right) . We are interested in that point on the critical surface that corresponds 
to the physical ratio of quark masses, i.e. R = mPj^^'^ /mFj^^^ — 24. From eq. 9 
we then find /ig/Tc ~ 1.25, which corresponds to /is ~ 'i.TaTc « 650 MeV and 
roughly coincides with the estimate given in ref. 12. Nevertheless our present 
uncertainties do not yet allow a reasonable error estimate on this value. 
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